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Abstract. It was recently shown that the theory of linear sto- 
chastic systems can be viewed as a particular case of the theory of 
linear systems on a certain commutative ring of power series in a 
countable number of variables. In the present work we study an 
interpolation problem in this setting. A key tool is the principle of 
permanence of algebraic identities. 
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1. Introduction 

There are numerous connections between classical interpolation prob- 
lems and optimal control and the theory of linear systems; see for 
instance [101 [I]- these settings, the coefficient space is the complex 
field C, or in the case of real systems, the real numbers M. Further- 
more, already from its inception, linear system theory was considered 
when the coefficient space is a general (commutative) field, or more 
generally a commutative ring; see p2l [25] . In [HJ [6] a new approach 



1991 Mathematics Subject Classification. 60H40, 93C05. 

Key words and phrases, white noise space, stochastic distributions, linear sys- 
tems on rings. 

D. Alpay thanks the Earl Katz family for endowing the chair which supported 
his research. 

1 



2 



D. ALPAY AND H. ATTIA 



to the theory of hnear stochastic systems was developed, in which the 
coefficient space is now a certain commutative ring 9^ (see Section |3] 
below). The results from |22l|25] do not seem to be directly applicable 
to this theory, and the specific properties of 5^ played a key role in the 
argumemts in [H [6] . 

The purpose of this work is to discuss the counterparts of classical 
interpolation problems in this new setting. To set the problems into 
perspective, we begin this introduction with a short discussion of the 
deterministic case. In the classical theory of linear systems, input- 
output relations of the form 

n 

(1.1) Un = ^hn--kUk, n = 0,l,... 

k=0 

where (M„)„eNo is called the input sequence, {yn)neNo is the output se- 
quence, and (/in)neNo is the impulse response, play an important role. 
The sequence (/in)nGNo consist of matrices (of common dimen- 
sions), and then the input and output sequences consist of vectors of 
appropriate dimensions. Similarly state space equations 

Un = CXn + DUn, n = 0, 1, . . . 

play an important role. Here x„ denotes the state at time n, and 
B, C and D are matrices with complex entries. The transfer function 
of the system is 

oo 

/i(A) = 5^/i„A", 

n=0 

in the case f ll.ip . and 

h{X) = D + XC{I-XA)-^B 

in the case of state space equations, when assuming the state at = 
to be equal to 0. Classical interpolation problems bear various ap- 
plications to the corresponding linear systems. See for instance (TUl 
Part VI], [21]. To fix ideas, we consider the case of bitangential inter- 
polation problem for matrix-valued functions analytic and contractive 
in the open unit disk (Schur functions), and will even consider only 
the Nevanlinna-Pick interpolation problem in the sequel to keep nota- 
tion simple, but it will be clear that the discussion extends to more 
general cases. Recall (see [TOl §18.5 p. 409]) that the bitangential in- 
terpolation problem may be defined in terms of a septuple of matrices 
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uj = C-, Att, A(^, B^, B^, r) by the conditions 

^ ResA=Ao(A/ - A^)-'B^S{X) = 
AoeB 

ResA=Ao5(A)C_(A/ - A^r' = C+, 
J2 ResA=Ao(A/ - A^)-'B+Si\)C4\I - A^y' = T, 

AqGD 

where A^ and have their spectra in the open unit disk, where 
(y4(^,i?+) is a full range pair (that is, controllable) and where (C_, A^) 
is a null kernel pair (that is, observable). We send the reader to [10] 
for the definitions. Moreover, F satisfies moreover the compatibility 
condition 

TA^ - Ac;T = B+C+ + B.C.. 
Let P be the matrix (see [IHl P- 458]) 



'1.2) P 



Pi T* 

r P2 



where Pi and P2 are the solutions of the Stein equations 

Pi - AiPiA^ = c*_c^ - 

P2 - A^PiAl = B+B\- B^B*,. 



Furthermore, and assuming the unknown function S to be C^^''- valued, 

J -- 



Ip 




When P is strictly positive, the solutions of the interpolation problem 
are given in terms of a linear fractional transformation based on a J- 
inner rational function G built from the septuple co via the formula (see 
[ini (18.5.6) p. 410]) 



e(A)^/+(A-Ao)(^^ il-XA,) 



(1.3) 

X P 



'C+ ~Bl\ f{XI -A^)-^ 

\i-\,Ai)-^ci -{I - x^Aiy^ci 

(,A^ - XoI)-'B+ {A^ - XoI)-'B_ 

where Aq is fixed on the unit circle and such that the various inverses 
exist in the above formula. An important fact is that the entries of Pi 
and P2 are rational functions of the entries of the matrices of u. As 
a consequence, there exists a rational function /(A), built from u and 
such that the entries of G are polynomials in A, with coefficients which 
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are themselves polynomials in the entries of the matrices of oj with co- 
efficients in Z. This fact will allow us in the sequel to use the principle 
of permanence of identities (see [9l p. 456]), to extend interpolation 
problem to a more general setting. 

Allowing in (11 .ip the input sequence {un)nmo consist of random vari- 
ables has been considered for a long time. On the other hand, allowing 
also the impulse response of the system to carry some randomness 
seems much more difficult to tackle. Recently a new approach to the 
theory of linear stochastic systems was developed using Hida's white 
noise space theory [TS], [12], [23], and Kondratiev's spaces of stochastic 
test functions and distributions [20]. In this approach, see [3], [5], [6], 
the complex numbers are replaced by random variables in the white 
noise space, or more generally, by stochastic distributions in the Kon- 
dratiev space, and the product of complex numbers is replaced by the 
Wick product. For instance, (II. ip now becomes 

n 

(1.4) yn = ^K-k()Uk, n = 0,l,... 

fc=0 

where the various quantities are in the white noise space, or more 
generally in the Kondratiev's space of stochastic distributions, and (> 
denotes the Wick product. An important role in this theory is played 
by a ring 9^ of power series in countably many variables with coefficients 
in C. This ring is endowed with a topology, which is that of the dual 
of a countably normed nuclear space. See Sections [2] and [31 Let us 
denote by 

(1.5) r(z) = ^v", 

an element of iH, where C denotes the set of sequences (ai,a2,...), 
whose entries are in 

No = {0,l,2,3,...}, 

and for which ^ Q for only a finite number of indices fc, and where 
we have used the multi-index notation 

The ring 91 has the following properties: 

(Pi) If r G $K and r(0, 0, 0, ■ ■ ■ ) 7^ 0, then r has an inverse in £H. 

(P2) If r e DV'''' is such that r(0, 0, 0, • • ■ ) = 0„xn and if / is a function 

of one complex variable, analytic in a neighborhood of the origin, with 



INTERPOLATION PROBLEM IN A COMMUTATIVE RING 



5 



Taylor expansion 



then, the series 



p=0 



oo 

dof. 



p=0 

converges in IH"^". Furthermore, if g is another function of one complex 
variable, analytic in a neighborhood of the origin, we have 

(1-6) (/^)(r) = /(r)^?(r). 

(P3) If r(^) = Eag£^a^" e 9^, then r*(^) =' E.g^^-^" ^ % where r* 
denotes the conjugate of the complex number r^. 

Property (Pi) implies in particular that a matrix A e j^^^x"^ is invertible 
in fK"^'" if and only if det A(0) ^ 0. This fact, together with (P2), 
allows to define expressions such as 

00 

(1.7) J^iX) = D + AC(/„ - AA)-^B = D + ^ A^CA'^-^B, 

k=l 

where A, B, C, and D are matrices of appropriate dimensions and 
with entries in SK, and where A is an independent complex variable. As 
explained in ^ this is the transfer function of some underlying linear 
systems, and is a rational function with coefficients in Di. 

The purpose of this paper is to explain how to tackle in the present 
setting counterparts of some classical interpolation problems which ap- 
pear in the theory of linear systems. To illustrate our strategy, we 
focus on the Nevanlinna-Pick interpolation problem, but our method 
works the same for the general bitangential interpolation problem. The 
computations done in the classical theory (that is, when the coefficient 
space consists of the complex numbers) extend to the case where C is 
replaced by the ring In some cases, such as Nevanlinna-Pick in- 
terpolation, this can be shown by direct computations. In the general 
case, one needs to use the principle of permanence of identities, see [HI 
p. 456]. We note that there are other commutative rings with proper- 
ties (Pi), (P2) and (P3) for which the above analysis is applicable. See 



The paper consists of six sections besides the present introduction. In 
the second section we review Hida's white noise space setting and the 
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Kondratiev spaces of stochastic distributions. The definition and main 
properties of the ring are given in Section [31 In Section H] we define 
and study analytic functions from an open set of C with values in Dl. 
In Section [5] we consider the Nevanlinna-Pick interpolation problem. 
In the last section we discuss the bitangential interpolation problem. 

2. The white noise space 

We here review Hida's white noise space theory and the associated 
spaces of stochastic distributions introduced by Kondratiev. See [IS], 
[in], [20], [23] • Let ^{R) denote the Schwartz space of smooth real- 
valued rapidly decreasing functions. It is a nuclear space, and by the 
Bochner-Minlos theorem (see p3| Theoreme 2, p. 342]), there exists a 

probability measure on the Borel sets B of the dual space 
such that 

(2.8) e ^ = / e*<"'^>rfP(w), Vs G ^(M), 

where the brackets (■, ■) denote the duality between ^(M) and ^'(M). 
The probability space 

W = (fi, B, dP) 

is called the white noise probability space. We will be interested in 
particular in L2(W), called the white noise space. For s G ^{M.), let 
Qs denote the random variable 

Qs{uj) = {u, s). 

It follows from (EH]) that 

||s||l2{k) = ||Qs||l2(w)- 

Therefore, Qs extends continuously to an isometry from L2(M) into 
L2(W). In the presentation of the Gelfand triple associated to the 
white noise space which we will use, we follow [20] ■ The white noise 
space L2(W) admits a special orthogonal basis {Ha)ae£, indexed by the 
set £ and built in terms of the Hermite functions hk and of the Hermite 
polynomials defined by 

CO 

H^{u) = l[h^,{Q~,^{u)). 

k=l 

We refer to f2U\ Definition 2.2.1 p. 19] for more information. In terms 
of this basis, any element of L2(W) can be written as 

(2.9) F = Y,faHo., UeC, 
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with 

There are quite a number of Gelfand triples associated to L2(W). In 
our previous works [2], [5], and in the present one, we focus on the one 
consisting of the Kondratiev space S\ of stochastic test functions, of 
W, and of the Kondratiev space S-\ of stochastic distributions. To 
define these spaces we first introduce for A; G N the Hilbert space Hfc 
which consists of series of the form (11. 5p such that 

(2.10) ||F|U =■ {y,{a^)^\f^\2nf'\ < oo, 

\a€i J 

and the Hilbert spaces "HJ, consisting of sequences {fa)aee such that 

iii^ir. =• fEi/«i'(2N)~'") ^ 

We note that, for F G "HJ, we have 

(2.11) lim ||F||; = 1/(0,0,0...)!', 

as can be seen, for instance, by applying the dominated convergence 
theorem to an appropriate discrete measure. Following the usage in the 
literature, we will also write the elements of lHy. as formal power series 
Sq6£ foMa- Note that {'Hk)keN forms a decreasing sequence of Hilbert 
spaces, with increasing norms, while ('H^)fcgN forms an increasing se- 
quence of Hilbert spaces, with decreasing norms. The spaces 5*1 and 
S^i are defined by the corresponding projective and inductive limits 

oo oo 

s^ = f^nk and = U ni 

k=l k=l 

The Wick product is defined on the basis {Ha)aei by 

It extends to an everywhere defined and continuous map from 5*1 x 5*1 
into itself and from S^i x S^i into itsel£| Let I > 0, and let > / + 1. 
Consider h GTi'i and u G "H^. Then, Vage's inequality holds: 



(2.12) \\hOu\\'k < A{k - l)\\h\ 



\u\ 



^The continuity properties are proved in [7, for a more general family of rings. 
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where 

(2.13) A{k -l)= f^(2N)('-'=)"^ < oo. 

See [201 Proposition 3.3.2 p. 118]. The following result is a direct 
consequence of (12.131) and will be useful in the sequel. 

Lemma 2.1. Let F eW^. Then, F^"" FO- - -OF G and 

n times 

(2.14) l|i^''^li;4-2<^(^(2)||F||;)", n = 1,2,3... 

Proof: We proceed by induction. The case n = 1 holds since 

< for Fen'^. 

Assume now that fl2.14p holds at rank n. Then, from f l2.12p we have 

I|fO("+i)||;^,<a(2)||f||;i|f^"II;^, 

<^(2))iiFii;^ (^(2)iiFii;)" 
^ (^(2)iiFii;)"^' 



A{2) 



3. The ring 9i 

The Kondratiev space S^i endowed with the Wick product is a com- 
mutative ring of sequences {ca)ae£, with properties (-Pi), (-P2) and (-P3), 
where in (Pi) one understands by evaluation at the origin the first co- 
efficient of the sequence. Using the Hermite transform (defined below), 
we view S-i as a ring of powers series in infinitely many variables. We 
point out that there are other commutative rings of sequences with 
properties (-Pi), (-P2) and (-P3), and for which a counterpart of the Her- 
mite transform holds. See [7]. 

The Hermite transform is defined by 

I{H^) = z", with aei and z = {z^, Z2, . . .) E C^. 

Then 
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It extends for F = J2aee^aHa € S-i by the formula I{F){z) = 
J2a£i^a^°'y converges in sets of the form 

Kp{R) = lzeC^ : J2 Izl^C^W <R^\, 

where p is such that F G Tip. The Kondratiev space S^i is closed under 
the Wick product, and we have 

IiFoG)iz) = I{F){z)I{G){z) and I{F + G){z) = I{F){z) + I{G){z) 

for any F,G & S-i. Therefore the image of the Kondratiev space S-i 
under the Hermite transform is a commutative ring, denoted by 

9^ =■ Im(l(S_i)). 

This ring was introduced in [6]. We transpose to it via the Hermite 
transform the properties of S^i. We have 

oo 

k=l 

We define the adjoint G* = (h,t) e O^'"^" of G = (gj,) G IH'"^"^ 
by hst{z) = gts{z) {t G {l,...,n} and s G {l,...,m}). Then for 
A G 9^"^™ and B G W''' we have 

(3.15) (AB)* = B*A*. 

Note that G*(0) = G(0)*, where G(0)* is the usual adjoint matrix. 

Definition 3.1. An element A G 9^"^" will he said strictly positive, 
A > 0, if it can he written as A = GG*, where G G ^^^^^ is invertihle. 
It will he said positive if G is not assumed to he invertihle. 

Lemma 3.2. Let A G 91"^". Then, A is strictly positive if and only 
if A{0) G C"^" is a strictly positive matrix (in the usual sense). 

Proof: If A = GG* with detG(O) ^ 0, then A(0) = G(0)G(0)* is 
a strictly positive matrix. Conversely, assume that A G 9^"^" is such 
that A(0) > 0. We write 

Aiz) = A(0) + (A(z) - A(0)) 

= vwnin + {^/my\M^) - a(o))(^/a(o))-^}a/a(o)- 

Let E{z) = ^A(0))-i(A(^)-A(0))(v/A(0))-i vanishes at z = (0, 0, . . .). 
Property (P2) with 

/(C) = 1 + - + + ■ ■ ■ = (1 + c)^/^ ICI < 1, 
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implies that A = C^, where C = /(E) = C*. □ 

Similarly, if A is positive, then A(0) is also positive, but the converse 
statement need not hold. Take for instance n = 1 and A{z) = Zi. Then 
it is readily seen that one cannot find r G such that Zi = r*{z)r{z). 

We define the ring of polynomials with coefficients in 9^ by 9^[A]. To 
avoid confusion between the variable A and the variables z we introduce 
the notation 

y{r) = r(0), re^R. 

Definition 3.3. A rational function with values in 9^"^™- is an expres- 
sion of the form 

(3.16) R(A)=p(A)(q(A))-i 

where p G (9^[A])"^™, and q G 9^[A] and such that J^(q(A)) ^ 0. 

Let R G 9^"^'"(A). Then, J^(R) G C"^'"(A), and it is readily seen that 

(3.17) (^(R))(A) = ^(R(A)). 

It is proved in [6] that every rational function with values matrices with 
entries in 91 and for which J^{c[{0)) ^ can be written as (11.71) . 

Example 3.4. Let r G £H. The function 

Fr(A) = (A - r)(l - Ar*)^^ G fH(A) 

is rational. It is defined for A G C such that 1 ^ A(j^(r))*. 

The next example of rational function need not be defined for A = 0. 

Example 3.5. Let r G 9^. The function 

Fr(A) = (A-r)(A-r*)-i G 9^(A) 

is rational. It is defined for A G C such that A ^ (^(r))*. 

4. Analytic functions with values in 9i 

It is possible to define analytic functions with values in a locally convex 
topological vector space (see for instance the discussion in [131 El [13 
[TB]). Here the structure of 9t allows us to focus, locally, on the classical 
definition of Hilbert space valued functions, as we now explain. 

Proposition 4.1. Let Q G C be an open set and let i : Q ^ be a 

continuous function. Then, f is locally Hilbert space valued, that is, for 
every (q G Q, there is a compact neighborhood K of Co (^nd a number 
Po such that i{K) C I^H' ). 
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Proof: Every (q & Q has a neighborhood K of the form Bs = {( & 
^ j ICo — CI ^ ^} for some 6 > 0. Since Bs is a compact set and f 
is continuous, f{Bs) is compact in and therefore strongly bounded. 
See [ni p. 54]. Thus there exists Po^N such that f(5^) G /(^pj and 
is bounded in the norm of liTi'p^^). See [121 Section 5.3 p. 45]. □ 

Therefore we can define an analytic function from i7 to 9^ as a continu- 
ous function which locally admits a power expansion with coefficients in 
one of the spaces I {'Hp). The following example shows that we cannot 
expect to have a fixed p in general. 

Example 4.2. Let i{X,z) = Then f is continuous (as a 

function of X) from C into fR, but there is no p such that f (A, z) (viewed 
now as a function of z) belongs to I ((Hp) for all A e C. 

Indeed, let Aq G C. We have 

oo oo 

(llf(Ao)li;)' = Yl l^'°l(2^)'' = 2-^5^n^^^°-^ < oo, 

n=l n=l 

for p > ReAo + 1. To show continuity at a point Aq G C, we take 
P > \Xo\ +2, and restrict A to be such that |Ao — A| < 1. Using the 
elementary estimate 

(4.1) |e^^ - e^^l < - Z2I ■ max |e^|, 

for Zi, Z2 G C, we have for n = 2, 3, . . . 

n2 — 2 < Inn e 2 

I I - 2 

and so 

00 

(l|f (A) - f (Ao) \U,)' = - I'-"'" 

n=2 

^ IA-AqP (lnri)2 

n=2 

and hence the continuity at the point Aq in the norm || • ||p+2, and hence 
in See in particular [12^ p. 57] for the latter. 

Recall that, in the case of Hilbert space, weak and strong analyticity are 
equivalent, and can be expressed in terms of power series expansions. 
The argument uses the uniform boundedness theorem. See Theo- 
rem VI. 4, p. 189]. We define the evaluation of an 9^- valued analytic 
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function at a point r G 9^. We first introduce 

where C C is open. 

Theorem 4.3. Let Q be an open subset ofC, and let i : Q be an 
analytic function. Let r G and let 

oo 

(4.2) f(c) = Ef"(c-^wr, 

n=0 

be the Taylor expansion around J^{r) G Q, where the f„ G "Hp^ /or some 
Po G N, anc? where the convergence is in Hp^- The series 

oo 

(4.3) f(r) = 5^f„(r-^(r))" 

n=0 

converges in H'^ for some q > po. 

Proof: Let i^' be a compact neighborhood of J^(r), and let po ^ N be 
such that i{K) C 'H'p^. Let furthermore R be the radius of convergence 
of the "Hpij-valued power series (14.21) . In view of (12.111) . there exists p, 
which we can assume strictly larger than po, such that 

(4.4) A(2)||r-^(r)||;<i?. 
On the other hand, using (12.141) . we obtain 

||f„(r - ^{r))X^, < A(2)||f„||;j|(r - J^{r))X^, 
<l|fn|i;o(^(2)l|r-^(r)||;)". 

In view of (14. 4p . the series 

oo 

Ellf'^llpo (^(2)l|r-^(r)||;)'^ 

n=0 

converges and so the series 



X]fn(r-^(r))^ 



n=0 



converges absolutely in /(Hp^g)- '-' 

The evaluation of f at r is defined to be f(r) given by (14. 3p . 

Proposition 4.4. We can rewrite the evaluation at r as a Cauchy 
integral 



2m J Q — Y 
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where the integration is along a circle centered at J^(r) and of radius 
r < R and in Q. 

Proof: As in Theorem 14.31 we consider a compact neighborhood K of 
J^{r), and let po be such that that i{K) C 'Hp^. We consider a simple 
closed path around J^(r) which lies inside K. 
We have 



dC 



2m J C-r 2m J C- .y (r) + J{y 

f(C) 




(C-^(r))"+i'^^ 

where we have used the estimates as in the proof of Theorem 14.31 and 
the dominated convergence theorem to justify the interchange of inte- 
gration and summation. □ 

Recall that a function / analytic and contractive in the open unit disk 
is called a Schur function. Furthermore, by the maximum modulus 
principle, / is in fact strictly contractive in D, unless it is identically 
equal to a unitary constant. We will call a function f analytic from 
the open unit disk D into 9^ a Schur function (notation: f G 5*9^) if the 
function 

A ^ ^(f(A)) 

is a Schur function. For instance, both 1 + ziz^ and 0.5 + lOzi — ?>z^ 
are Schur functions. We now define the analog of the open unit disk 

by 

= {r G fH; =^(r) G D}, 
and the analog of strictly contractive Schur functions as the set of 
analytic functions f : D — )• such that the function A c^(f(A)) is a 
strictly contractive Schur function. 

Theorem 4.5. i E S<}\ is a strictly contractive Schur function if and 
only if { : D — 7- Dlo is analytic. 

Proof: If f is analytic from D into and such that the A t-)- J^(f(A)) 
is a strictly contractive Schur function, it means by definition that the 
range of f lies inside ^Hb- Conversely, let f : D — )■ be analytic and 
such that J^{i) is a strictly contractive Schur function. Then for every 
< r < 1, there exists /c G N (which may depend on r) such that 
f(|A| < r) C /("H'fe). We can write f as 
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f(A) = ^A'^f„, 

n=0 

where |A| < r and fn G /(H^,). Now ^(f)(A) = Er=o^"^(fn) for 
|A| < r. Since this holds for all r e (0, 1) the function A t— ?■ f(A) has 
range inside D^b- D 



5. Nevanlinna-Pick Interpolation 

In this section we solve the following interpolation problem {IP)- 

Problem 5.1. Given ri G N and points ai, . . . , a„, bi, . . . , b^i G Din, 
find all Schur functions f with coefficients in Dl such that f(ai) = bj 
for i = 1,2, . . . ,n. 

The solution of this problem under the assumption that some matrix 
is strictly positive, is presented in Theorem 15.31 below. We first give 
some preliminary arguments, and note that if f is a solution of the 
interpolation problem 15. then / = J^(f) is a solution of the classical 
interpolation problem 

(5.5) f{ai) = bi, i= l,...,n, 

where we have set 

(a-j) and bi = J^(bj), i = 1, . . .n. 



This last problem is solved as follows: let P denote the nxn Hermitian 
matrix with ij entry equal to 

1 - bib* 

(5.6) 

A necessary and sufficient condition for (15. 5p to have a solution in the 
family of Schur functions is that P > 0. We will assume P > 0. Set, 
in the notation of the introduction, 

Al = A = diag {al,a*2 . . . ,<), 

- f^A = ^ ■ ■ ■ ^ ^ '^"^ 

(5.7) \B_ -\b\ 6* ■■• bl 



J 



1 
-1 



Furthermore, specializing the formula for B given in the introduction 
with 2;o = 1, or using the formula arising from the theory of reproducing 
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kernel Hilbert spaces (see [II],[T]), set 

e(A) = h-{i- x)c{i^ - XA)-'p-\i - Arc* J ^1^]) . 

We now gather the main properties of the matrix-valued function G 
relevant to the present work. For proofs, we refer to [1], [10], jllj . 

Proposition 5.2. The following hold: 

(a) The function G is J -inner with respect to the open unit disk. 

(b) G has no poles in D and a{X)a + 6(A) 7^ for a// A G D and all a 
in the closed unit disk. 

(c) The identity 

(5.8) (1 -6,)G(a,)=0, i = l,...n. 

is valid, (d) The linear fractional transformation 

i i^-,^ dcf. a(A)or(A) + &(A) 
= c(A)a(A) + rf(A) 

describes the set of all solutions of the problem (15. 5p in the family of 
Schur functions when o varies in the family of Schur functions. 



To solve the interpolation problem 15. II we introduce the matrices A, C 
and P, with entries in 9^, built by formulas (15. 6 p and (15. 7p . but with 
ai, . . . , bi, . . . , b„ instead of Oi, . . . , 61, . . . , 6„. Note that P > 
since P > 0, and we can define the 9^^^ ^-valued function as G but 
with A, C and P instead of A, C and P. We have 

,y{A) = A, y{C) = C, and ^(P) = P. 

Furthermore, 

(5.9) ^(e(A)) = G(A). 

Theorem 5.3. Assume P > 0. Then, there is a one-to-one correspon- 
dence between the solutions f of the problem \5.1\ in S<}\ and the elements 
g G via the linear fractional transformation f = T©(g). 

Proof: We first claim that the matrix-valued function satisfies the 
counterparts of (15. 8p . that is, 

(5.10) (1 -b,)e(ai) = 0, i = l,...n. 

This is done using the permanence of algebraic identities. See ^ p. 
456] for the latter. Indeed, the function 

(det(J„-AA))(det(J„-A*))(detP)( J] (1 - a,a*))G(A) 
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is a polynomial in A with coefficients which are themselves polynomials 
in the Oj and the bj, with entire coefficients. Therefore, multiplying 
both sides of (15. 8p by the polynomial function 

(det(4-AA))(det(/„-A*))(detP)( J] (1 - a,a*)) 

e,j=l,...n 

evaluated at A = Oj (i = 1, 2, . . . , n), and taking the real and imaginary 
part of the equalities (15. 8p . we obtain for each i four polynomial identi- 
ties in the 4n real variables Rea^, Kebj, Ima^, Imbj, with j = 1, . . .n, 
with entire coefficients, namely 

Re I (1 -bi) det (/ - a,A) det (/ - A*) det ^ JJ ~ aea*)Q{ai) I 

I e,j=l,...n J 



(0 0), 

- a£a*)Q{ai 



Im I (1 -bi) det (/ - aiA) det (/ - A*) det ^ JJ " 
I e,j=i,-n 

= (0 0). 

It follows (see p. 456]) that these identities hold in any commutative 
rings, and in particular in Di: 

Re I (1 -bi) det (/ - a^A) det (/ - A*) det P JJ ~ a^a*)0(ai) I 
I e,j=i,-n J 

= (0 0), 

Im I (1 -hi) det (/ - a^A) det (/ - A*) det P JJ ~ a^a*)0(ai) I 

I e,j=l,...n J 

= (0 0). 

We now use the fact that we are in the ring 9^. Because of the choice 
of the SLj, the element 

det (J - a,A) det (J - A*) JJ (1 - a^a*) 

e,j=l,...n 

is invertible in £H. When furthermore P > we can divide both sides 
of the above equalities by 

det (/ - a^A) det (/ - A*) JJ (l-a^a*)detP 

e,j=l,...n 

and obtain (15.101) . 
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Let now r G SfR, and let u, v be analytic valued functions defined by 



Using flS.lOp we have that 

u(aj) = biv(ai), i = l,...,n. 

To conclude, we need to show that v(aj) is invertible in £H for i = 
1, . . . n. But we have 

J^(v(ai)) = c{ai)y{r){ai) + d^ai), i = l,...,n. 

Since the function 6(A) = J^(©(A)) is J-unitary on the unit circle 
and has no poles there. Therefore, we have c{ai)J^{r){ai) + d{ai) ^ 
(see item (6) in Proposition 15.21) . and hence v(aj) is invertible in $H. 
Therefore uv~^ = T©(r) is a solution of the interpolation problem. 

Assume now that f is a solution. Then, we know from the discussion 
before the theorem that there exists a Schur function cr(A) such that 



Comparing with (15.111) . we obtain J^(r(A)) = cr(A), and hence r G Syi. 

□ 

6. More general interpolation problem 

The matrix-valued function defined by (II. 3p and describing the set 
of solutions of the bitangential problem satisfies the conditions 



See also [1]. As for the Nevanlinna-Pick case, these conditions can 
be translated into a finite number of polynomial equations with co- 
efficients in Z, and the principle of permanence of identities allows to 
extend these properties in the case of a commutative ring. On the other 
hand, we do not know how to extend the third interpolation property. 





^ ResA=Ao ((A/ - 5-)e(A) = 



AqGD 
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and so the method is not apphcable to the most general bitangential 
interpolation problem. On the other hand, if we restrict the parameter 
to be a constant contractive matrix, the third condition also translates 
into a polynomial identity with entire coefficents, and the same method 
can still be used. The case of functions with poles inside the open unit 
disk, or the degenerate cases, are more difficult to treat, and will be 
considered elsewhere. 
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